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We discuss the properties of a self-organized critical forest- 
fire model which has been introduced recently. We derive 
scaling laws and define critical exponents. The values of these 
critical exponents are determined by computer simulations in 
1 to 8 dimensions. The simulations suggest a critical dimen- 
sion d c = 6 above which the critical exponents assume their 
mean-field values. Changing the lattice symmetry and allow- 
ing trees to be immune against fire, we show that the critical 
exponents are universal. 
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I. INTRODUCTION 

Some years ago, Bak, Tang, and Wiesenfeld intro- 
duced the sandpile model which evolves into a critical 
state irrespective of initial conditions and without fine 
tuning of parameters [1]. Such systems are called self- 
organized critical (SOC) and exhibit power-law correla- 
tions in space and time. The concept of SOC has at- 
tracted much interest since it might explain the origin of 
fractal structures and 1/ f -noise. Models for earthquakes 
[2,3], the evolution of populations [4,5], the formation of 
clouds [6] and river networks [7], and many more have 
been introduced and investigated by computer simula- 
tions, but in contrast to the sandpile model, most of 
these SOC models are barely understood. For most of 
them exists no proof that they really are critical, and lit- 
tle analytic treatment has been presented so far. It even 
has been conjectured that systems without conservation 
laws cannot be SOC. 

This conjecture has been refuted when a forest-fire 
model without conservation laws was shown to be crit- 
ical under the condition that time scales are separated 
[8]. In one dimension, where the model is still nontrivial, 
the exact values of the critical exponents have been cal- 
culated [9], thus proving the criticality of the model. A 
scaling theory presented in [8] suggested classical values 
for the critical exponents in two and more dimensions but 
has been shown to be too simple [10-12]. Scaling laws 
for the static properties of the model have been given 
in [12]. The values of several critical exponents for the 
two-dimensional model have been determined by com- 
puter simulations [10,12]. These simulations differ in the 
values of two of the critical exponents. Simulations in 
up to 6 dimensions [11] suggest that the critical behavior 



of the forest-fire model becomes percolation-like above 
6 dimensions. 

In this paper, we present a scaling theory for the statics 
and dynamics of the system and present simulation re- 
sults in 1 to 8 dimensions. Additional critical exponents 
and scaling relations are introduced, and the values of 
the critical exponents are determined by computer sim- 
ulations, which are performed closer to the critical point 
than the earlier simulations. The universality of the val- 
ues of the critical exponents is investigated by changing 
the lattice symmetry and by considering the case of non- 
vanishing immunity. 

The outline of the paper is as follows: In Sec. II, the 
rules of the model are introduced and the origin of the 
SOC behavior is explained. In Sec. Ill, the scaling theory 
of the model is presented. In Sec. IV, simulation results 
in 1 to 8 dimensions are given, and the critical behavior 
in high dimensions is discussed. Sec. V investigates the 
universality of the critical exponents. In the final section, 
the results are summarized and discussed. 



II. THE MODEL 

The forest-fire model is a stochastic cellular automaton 
which is defined on a rf-dimensional hypercubic lattice 
with L d sites. Each site is occupied by a tree, a burning 
tree, or it is empty. During one time step, the system is 
parallely updated according to the following rules 

• burning tree > empty site 

• tree > burning tree, if at least one nearest neigh- 
bor is burning 

• tree — > burning tree with probability /, if no 
neighbor is burning 

• empty site — > tree with probability p. 

An even more general forest-fire model also contains an 
immunity g [13,14]. In its original version, introduced by 
P. Bak, K. Chen, and C. Tang, the forest-fire model con- 
tained only the tree growth parameter p [15]. This ver- 
sion of the model shows regular spiral-shaped fire fronts 
in the limit of slow tree growth [16,17]. Throughout this 
paper, we will assume that the system size L is large 
enough such that no finite-size effects occur. In the sim- 
ulations, we have always chosen periodic boundary con- 
ditions. 
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Starting with arbitrary initial conditions, the system 
approaches after a transition period a steady state the 
properties of which depend only on the parameter values. 
Let p e , p t , and pj be the mean density of empty sites, 
of trees, and of burning trees in the steady state. These 
densities are related by the equations 

Pe + Pt + PS = 1 (2-1) 

and 

Pj=PPe- (2.2) 

The second relation says that the mean number of grow- 
ing trees equals the mean number of burning trees in the 
steady state. Large-scale structures and therefore criti- 
cality can only occur when the fire density pj is small. 
When the fire density is large, trees cannot live long 
enough to become part of large forest clusters. Eq. (2.2) 
shows that the fire density becomes small when the tree 
growth rate p approaches zero. In addition, the lightning 
probability / must satisfy 

f<p. (2.3) 

Otherwise a tree is destroyed by lightning before its 
neighbors grow, and no large-scale structures can be 
formed. 

These conditions are not yet sufficient to bring about 
critical behavior in the forest-fire model. When light- 
ning strikes a small forest cluster, it burns down very 
fast, before any tree can grow at its edge. But when 
lightning strikes a large forest cluster, it needs some time 
to burn down, and new trees might grow at the edge of 
this cluster while it is still burning so that the fire is never 
extinguished. In order to observe critical, i.e. self-similar 
behavior, small and large forest clusters must burn down 
in the same way. We therefore choose the tree growth 
rate p so small, that even the largest forest cluster burns 
down rapidly, before new trees grow at its edge. In this 
case, the dynamics of the system depend only on the ra- 
tio f/p, but not on / and p separately. When / and p are 
both decreased by the same factor, the overall time scale 
of the system is also changed by this factor, but not the 
number of trees that grow between two lightnings and 
therefore not the size distribution of forest clusters and 
of fires. The condition that forest clusters burn down 
rapidly can be written in the form 

p<T- 1 (s max ), (2.4) 

where T(s max ) is the time the fire needs to burn down 
a large forest cluster and will be determined below (see 
Eq. (3.20)). In the simulations, this condition is most eas- 
ily realized by assuming that forest clusters burn down 
instantaneously, i.e. during one time step. This simpli- 
fied version of the SOC forest-fire model [8] has been 
mentioned first in [18]. 

The inequalities (2.3) and (2.4) represent a double sep- 
aration of time scales 



T(s max ) <. p- 1 <. f- 1 , (2.5) 

which is the condition for SOC behavior in the forest-fire 
model. The time in which a forest cluster burns down 
is much shorter than the time in which a tree grows, 
which again is much shorter than the time between two 
lightning occurrences. Separation of time scales is quite 
frequent in nature, while the tuning of parameters to a 
certain finite value only takes place accidentally. Thus, 
the forest-fire model is critical over a wide range of pa- 
rameter values. A snapshot of the critical state is shown 
in Fig. 1. 




FIG. 1. Snapshot of the SOC state in 2 dimensions. 
Trees are black, empty sites are white. The parameters are 
L = 1024 and f/p = 1/500. 

III. SCALING LAWS AND CRITICAL 
EXPONENTS 

In this section, we will derive scaling laws and relations 
between the critical exponents for the SOC forest-fire 
model. 

First, we calculate the mean number s of trees that are 
destroyed by a lightning stroke. During one time step, 
there are 

fptL d 

lightning strokes in the system and 

PPeL d 

growing trees. In the steady state, the number of growing 
trees equals the number of burning trees, and therefore 
the mean number of trees destroyed by a lightning stroke 
is 

_ pp e p 1 - p t , , 

s= — ~- . (3.1) 

fPt f Pt 
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In the last step, we have neglected the fire density which 
is very small due to time scale separation. For small 
values of f/p, the forest density p t assumes a constant 
value. If this constant value is less than 1, the second 
factor on the right-hand side of Eq. (3.1) is also constant 
for small f/p, and Eq. (3.1) then represents a power law 

stxtf/p)- 1 . (3.2) 

In d > 2 dimensions, the critical forest density 

Pt = I™ Pt, (3-3) 

in fact, must be less than 1, as the following considera- 
tion indicates: If the critical forest density were p c t = 1 
in d > 2 dimensions, p t would be very close to 1 for small 
values of f /p. Then the largest forest cluster would con- 
tain a nonvanishing percentage of all trees in the system, 
and the average number s of trees burned by a lightning 
stroke would diverge in the limit L — ► oo with fixed f/p, 
in contradiction to Eq. (3.1). In one dimension, there is 
no infinitely large forest cluster in the system as long as 
Pt < 1, and therefore the critical forest density is p'\ = 1. 
Nevertheless Eq. (3.2) holds also in 1 dimension since 
the forest density approaches its critical value only log- 
arithmically slowly, as will be shown below. Eq. (3.2) 
indicates a critical point in the limit f/p — ► 0. Close 
to this critical point, i.e. if / <C P, there is scaling over 
many orders of magnitude. 

Let n(s) be the mean number of forest clusters per 
unit volume consisting of s trees. Then the mean forest 
density is 



A = l/(3-r). 



(3.8) 



Pt 



\sn(s), 



(3.4) 



and the mean number of trees destroyed by a lightning 
stroke is 



J2 s ' 2n ( s )/Pt 



(3.5) 



In the case r = 2, the right-hand side of Eq. (3.6) 
acquires a factor l/ln(s max ) and reads now 



n(s) cx s T C(s/s max )/ln(s max ) 



(3.9) 



since the forest density given by Eq. (3.4) must not di- 
verge in the limit f/p — ► 0. The mean number of forest 
clusters per unit volume Ei° n(s), therefore, decreases 
to zero for f/p — ► 0, and consequently the forest density 
approaches the value 1. 

We also introduce the cluster radius R(s) (radius of 
gyration) which is the mean distance of the trees in a 
cluster from their center of mass. It is related to the 
cluster size s by 



s cx R(sf 



(3.10) 



with the fractal dimension p. 

The pair connectedness C(x, f/p) is the probability 
that a site at distance x from an occupied site is also 
occupied and belongs to the same cluster [19]. The cor- 
relation length £ is defined by 



E x x 2 C(x,//p) 

ExC(x,// P ) 

Er s »( s ) E*=i 1 

_ J2T 2sn( s ) -sR 2 (s) 

Er sn ( s ) ■ s 

/oo 
dss 2- T+ 2/, C ( s/Smax) 

oc (//P)" 2A/ ^, (3.11) 

where (. . .) denotes the average over all clusters of size s. 
We conclude 

£ cx (f/p)~ v with v = X/p. (3.12) 

Another quantity of interest is the mean cluster radius 



Since lim^ / p ^o pt is finite and s diverges cx (f/p) 1 , these 
equations imply that n(s) decreases at least like s~ 2 but 
not faster than s~ 3 . As long as the system is not exactly 
at the critical point f/p = 0, i.e. for nonvanishing f/p, 
there must be a cutoff in the cluster size distribution for 
very large forest clusters. We conclude that [8] 



n(s) cx s T C(s/s max ) 
with 2 < r < 3 and 



(f/p) cx (f/p) 



" A ocs A . 



(3.6) 



(3.7) 



The cutoff function C(x) is essentially constant for x < 1 
and decreases to zero for large x. Eqs. (3.5) - (3.7) yield 
s cx s^~^, which leads to the scaling relation 



R = sn(s)R(s)/Y^ sn(s) 

s=l s=l 

' cx (f/p)- ( "- x+1) , if v - X + 1 > ; 
= const., ii v — A + 1<0. 

This leads to 

R cx (f/p)-" with v = max(0, v - (X - 1)) . (3.13) 

In percolation theory, the hyperscaling relation 

d = p(r-l) (3.14) 

is satisfied, but it is not satisfied in the SOC forest-fire 
model in d = 2, as first stated in [12], where also an 
interpretation of this relation is given: If Eq. (3.14) is 
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satisfied, every box of l d 1 sites contains a spanning 
piece of a large cluster when the system is at the critical 
point. In the forest-fire model, there are at least in d = 2 
many regions which contain no large forest cluster (see 
Fig. 1), and consequently d < p(r — 1). 

The mean forest density p t approaches its critical value 

= liirif/p^o P~t v i a a P ower law 

P e t -pt<x(f/p) 1/e . (3.15) 
One might try to calculate 1/6 in the following way 

CO 

p r t - p t cx ^2 s 1 ~ T (1 - C(s/s max )) 

s = l 

dss 1 ~ T (1 - C(s/s max )) 

«^-i«(//p)( T - 2 ^ 3 - T ). (3.16) 
The simulation results for r and 6, however, show that 

(r-2)/(3-r) 
is much smaller than 1/8 (see Tab. I). We conclude that 



dxx 1 ~ T (C(x) - 1) = 



(3.17) 



and that the above formula Eq. (3.6) for n(s) acquires an 
additional factor 

1 - (f/p) 1 ' 6 An(s) . 

Eq. (3.17) says that the number of trees in the system 
is not influenced by the cutoff" function. All trees that 
were in large forest clusters if there were no cutoff can be 
found in smaller clusters. The size distribution of forest 
clusters therefore has a bump (see Fig. 2). This expla- 
nation of the bump has already been suggested in [10]. 
The function An(s) determines how the distribution of 
forest clusters of size s < s max deviates from the critical 
distribution. Since the mean density of empty sites de- 
creases with decreasing f/p, less small forest clusters are 
produced when the system is closer to the critical point 
f/p — ► 0. Consequently An(s) is negative for small s. 
Since the forest density increases with decreasing f/p, 
the function An(s) must be positive for larger values of 
s, leading to 



dss T C(s/s max ) An(s) > . 



Finally, we introduce dynamical exponents character- 
izing the temporal behavior of the fire. Let T(s) be the 
average time a cluster of size s needs to burn down when 
ignited, and N(T) the portion of fires that live exactly for 
T time steps. Then the exponents b and p' are defined 

by 



From 

N{T)dT cx sn(s)ds 
follows the scaling relation 

b = p'(T -2) + l. (3.19) 
The time scale of the system is set by 

7m ax = T(s max ) oc (f/p)-"' with v' = X/p'. (3.20) 
The dynamical critical exponent z is defined by 

which leads with (3.12) and (3.20) to 

z = v'/v = p/p'. (3.21) 

The condition of time scale separation now can be ex- 
pressed in terms of the critical exponents and reads 



(f/p)-"' <p-' <f~\ 
or equivalently 

/ <p<r' /(1+i/) . 

The average lifetime of fires is 



(3.22) 



(3.23) 



f = J2 sn(s)T(s)/ E sn ( s ) K (//P)" 5 ' ( 3 - 24 ) 



with 



v' = max(0,i/ - (A - 1)) . 



(3.25) 



The average number N s (t) of trees that burn t time steps 
after a cluster of size s is struck by lightning enters the 
definition of the temporal fire-fire correlation function 
G(r) 



CO CO 



G(T)(xj2 n ( s ) s J2 Ns ^ Ns ( t + TS >- ( 3 - 26 ) 



s=l t=0 

The power spectrum is the Fourier transform of the fire- 
fire correlation function 

G(ui) = 2 / dr G(t) cos(wt) cx ui~ a for small ui. 
Jo 

(3.27) 

Using the function N s (t), the mean chemical distance of 
the trees in a cluster from the site of the lightning stroke 
can be defined [20]: 



s cx T(sf and N(T) cx T~ 



(3.18) 



(3.28) 
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This definition is analogous to the definition of the clus- 
ter radius as the mean distance of the trees of a cluster 
from the center of mass. When the structure of the for- 
est clusters is compact or fractal, the chemical radius is 
proportional to the life time of the fire: 

t(s) oc T(s) , (3.29) 

and i(s) introduces no new critical exponents. Conse- 
quently T max is proportional to a chemical correlation 
length, which could be defined by a relation similar to Eq. 
(3.11). Our computer simulations confirm Eq. (3.29). 

IV. VALUES OF THE CRITICAL EXPONENTS 
IN 1 TO 8 DIMENSIONS 

In this section, we determine the values of the critical 
exponents in 1 to 8 dimensions. In one dimension, the 
critical exponents can be determined analytically, as was 
done in [9]. In higher dimensions one has to resort to 
computer simulations. Here, we shortcut the exact eval- 
uation of the critical exponents in d = 1 by using simple 
arguments. In subsection IV B we present the simulation 
results in 2 and more dimensions [21]. In subsection IV C 
we discuss the behavior of the system on a Bethe lattice. 

A. The critical exponents in one dimension 

In one dimension, the critical forest density p'\ equals 
1, since otherwise there were no infinitely large forest 
cluster in the system. The consideration after Eq. (3.8) 
shows that consequently r = 2 and (via scaling relation 
Eq. (3.8)) A = 1. In the steady state, the density of forest 
clusters ^2n(s) is constant, and therefore [22] 

oo 

J2n( S ) = (l-p t -(f/p)p t )/2, 

s = l 

which leads together with Eq. (3.9) to 

(1 - p t ) oc l/ln(s max ) 

and 1/8 = 0. One-dimensional forest clusters are com- 
pact, therefore p = 1 and (with Eq. (3.12)) v = 1. From 
T(s) oc R(s) it follows p' = v' = z = 1. The exact calcu- 
lation in [9] yields additional logarithmic corrections: 

oc £ oc T max oc (p//)/ln(p//). 

The fourier transform of the temporal correlation func- 
tion is [9] 

G(lo) oc lo~ 2 (1 + const. • ln(ws max )) , (4-1) 

indicating a deviation from the trivial w _2 -dependence 
towards 1/w-noise. Tab. I summarizes the values of the 
critical exponents. They are confirmed by our simula- 
tions. 



B. Simulation results in 2 and more dimensions 

We obtained the values of the critical exponents in d > 
2 dimensions by computer simulations using the same 
method as in [10] which iterates the following rules: 

1. Choose an arbitrary site in the system. If it is not 
occupied by a tree, proceed with rule 2. If it is 
occupied by a tree, then ignite the tree and burn 
down the forest cluster to which the tree belongs. 
While burning the trees, evaluate the size and the 
radius of the cluster, the lifetime of the fire and the 
function N s (t). Proceed with rule 2. 

2. Choose p/ f arbitrary sites in the system and grow 
a tree at all chosen empty sites. Proceed with rule 
1. 

Time scale separation is perfectly realized by these rules, 
and Eq. (3.1) is satisfied. 

In order to assure that the system is in the steady state, 
a sufficiently large number of time steps have been dis- 
carded in the beginning of each simulation. The simula- 
tions were performed on a DECstation 5000. The system 
length L for each dimension used in our simulations is 
given in Tab. I. The simulations could not be performed 
arbitrarily close to the critical point f /p — ► 0, since the 
correlation length £ diverges, and, therefore, finite-size 
effects occur when £ exceeds the system size. Especially 
the distance-related exponents p, v and v could not be 
determined with available computer capacity in high di- 
mensions. In the following, we discuss the simulation re- 
sults obtained for the critical exponents and the critical 
forest density. 

a. Exponent r In contrast to the rest of the criti- 
cal exponents, the exponent r can also be determined in 
higher dimensions, since the cluster distribution "feels" 
more the overall system size L d than the edge length L. 
Therefore the exponent r could be measured in 1 to 8, 
12 (L = 4) and 16 (L = 3) dimensions (see Fig. 2). Since 
we evaluated the forest clusters struck by lightning, we 
measured the fire distribution sn(s), which gives r — 1. 
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FIG. 2. Normalized cluster distribution n(s)/n(\) in d = 2 
to 6 dimensions. The values of f/p are 1/32000, 1/2000, 
1/1000, 1/250 and 1/125 from right to left. The negative 
slope yields the exponent r. 

The values of f/p which just don't lead to finite-size 
effects are determined by means of the integrated fire dis- 
tribution P(s), which is defined [10] by 

/>oo 

P(s) = / ds's'n(s'). (4.2) 

J s 

A plot of P(s)/P(l) for different values of f/p in d = 6 
is shown in Fig. 3. 




10° 10' ~~ 10 4 10' 

s 

FIG. 3. Normalized integrated fire distribution P(s)/P(l) 
for f/p = 2/125, 1/125, 1/250, 1/500 in d = 6 (from left to 
right). 

As long as the finite size of the system is not impor- 
tant, the scaling regime of P(s) becomes larger, if one 
approaches the critical point. As soon as P(s) feels the 
finite system size, the curves bend upwards (see the two 
uppermost curves in Fig. 3). Evaluating the curves not 
affected by finite-size effects, we obtain the values of r 
shown in Tab. I. For d = 2, they are consistent with 
[10,12]. The values of r for d = 2 to 6 are compatible 
with those given in [11] but have a smaller error since 
they have been obtained using larger systems closer to 
the critical point. Above d = 6 dimensions, the numeri- 
cal results suggest r = 2.5 as in percolation theory. Our 
results indicate that the SOC forest-fire model has the 
upper critical dimension d c = 6, as already conjectured 
in [10,11]. An analytic proof, however, is still missing. 

b. Exponent A The determination of A is slightly 
more difficult. Consider again the normalized integrated 
distribution function P(s)/ P(l): 

j p(s) = } r^/- r c( s y Smai ) 

n p(1) - 11 f™ d8 > 8 > 1 - T C( 8 >/ 8max ) 

{• CO fOQ 

= In / dx x l ~ T C(x) — In / dxx l ~ T C(x) 

^s/s max -/l/s mas 

= -4(s/s m ax) - A(l/s max ). (4.3) 



Eq. (4.3) implies that the plots of ln(P(s)/P(l)) vs. ln(s) 
for different values of f/p can be made to coincide, when 
they are shifted horizontally and vertically (see Fig. 4). 




FIG. 4. Normalized integrated fire distribution P(s)/P(l) 
for f/p = 1/16000 and f/p = 1/32000 in d = 2. 

The logarithm of the horizontal scaling factor divided 
by m ((//p)i/(//p)2) then yields A. Our simualtion re- 
sults in d = 2 and 3 are given in Tab. I. They satisfy the 
scaling relation Eq. (3.8). 

With increasing dimension, the range of admissible val- 
ues of f/p decreases due to finite-size effects. Conse- 
quently, the error in A becomes fairly large in d > 4. 
Hence, in these dimensions, its value was calculated us- 
ing the scaling relation Eq. (3.8). 

For d = 2, the value of A has also been determined 
in [12] which agrees with ours, and in [10] which does 
not agree with ours and is inconsistent with the scaling 
relation Eq. (3.8). 

c. The critical forest density p'\ and the exponent 8 
In our simulations the forest density p t was determined 
by dividing the number of burnt-down clusters by the 
number of iterations. In order to calculate the critical 
forest density p'\ and the exponent 8, one needs at least 
three simulations at different f/p. In d = 2,3 we fitted 
the results pt(f /p) with the function pi — const. • (f/p) 1 / 6 
(Fig. 5). 
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FIG. 5. pi - p t in d = 2 (o) and d = 3 (□). The inverse 
slope yields the exponent 8. 

The error in 1/8 is small enough to rule out the pos- 
sibility that 1/8 = (t - 2)/(3 - r) (see Eq. (3.16)). In 
d > 3 the change in pt(f /p) and the error in p t are of the 
same order of magnitude, making the fit useless. 

The results for p'\ and 1/8 as well as 1/8 in perco- 
lation theory and (r — 2)/(3 — r) are given in Tab. I. 
Also shown in Tab. I are the corresponding values of the 
mean-field forest density. In mean-field theory the crit- 
ical forest density equals l/(2d — 1), when taking into 
account that burning trees always have an empty neigh- 
bor where the fire has come from. p'\ is indistinguishable 
from the mean-field density for d > 5. 

Our value of p'\ in d = 2 is consistent with [10,11]. 
The result in [12] for d = 2 as well as the results in [11] 
for d = 3 to 6 are moderately larger than ours, which 
is probably due to the smaller system size of these other 
simulations. 

Our values of 1/8 are compatible with [10,11], whereas 
the result in [12] for d = 2 seems to be too small. In 
three dimensions our result for 8 is already very close to 
the percolation value. 

d. Exponents p, v and v The fractal dimension p 
of the forest clusters is obtained from the slope of the 
cluster radius R(s) (see Fig. 6). 
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FIG. 6. Cluster radius R(s) in 2 to 6 dimensions (from left 
to right). The inverse slope yields the fractal dimension fi. 

The distance between two lattice sites x, x' is defined 

by 

d 

d(x,x') = |, (4.4) 

8 = 1 

since this is the length of the shortest path between them. 
In 6 and higher dimensions the i?(s)-plot has no scaling 
region, therefore p could only be determined in d = 2 to 
4 and, with relatively large error, in d = 5. Unlike [8,10], 
we find p < d in d = 2, in accordance with [12]. From the 
results in Tab. I it seems that p approaches the fractal 



dimension of percolation clusters and we suggest that 
H = p peTC for d > 6. The hyperscaling relation Eq. (3.14) 
is definitely violated in d = 2, but holds within error 
margins in d = 3, 4, 5. 

The correlation length is dominated by large clusters 
and consequently large radii. Therefore the exponent v 
could directly be determined only in d = 2 dimensions. 
Since R is less divergent than £, the exponent v could 
also be determined in d = 3 dimensions. Our value of v 
in d = 2 is compatible with [10,12], and all results (see 
Tab. I) satisfy the scaling relations Eqs. (3.12) and (3.13). 

e. Exponents p' , v' , and v' The chemical dimension 
p' [20] can be determined even in d = 16, since the scaling 
regime extends over at least two decades. In 3 and more 
dimensions, the results are compatible with the mean- 
field value p' = 2. Just as v and v, the exponents v' 
and v' could only be determined in d = 2 and d = 2, 3 
respectively. The results (see Tab. I) obey the scaling 
relations Eqs. (3.25) and (3.20). 

/. Exponents b and z Knowing the exponents p, p' 
and r, one can easily calculate the exponent of the distri- 
bution of fire lifetimes b (see Eq. (3.19)) and the dynami- 
cal exponent z, which characterizes the relation between 
the time scale and the length scale of the system (see 
Eq. (3.21)). The results are shown in Tab. I. Both expo- 
nents seem to approach the value 2 for d — ► 6. 

g. Exponent a Assuming a simplified law N s (t) oc 

-i 0(J 1 ( S ) — 1) ; one can show that (to leading order in 
the frequency) a = 2 for d > 6. For 1 < d < 6, a has 
to be determined by simulations. The results are shown 
in Tab. I. The non-monotonous behaviour of a between 
d = 1 and d = 6 is not so surprising if one recalls the 
intricate dependence of a on p' and r (see Eq. (3.26)). 

In summary, we can conclude that the SOC forest- 
fire model is likely to have an upper critical dimension 
d c = 6, above which the critical exponents are identical 
with those of mean-field-theory, which again is identical 
to the mean-field-theory of percolation. The strongest 
evidence for this behaviour comes from the exponent r, 
which approaches the percolation value r perc = 5/2 for 
d — ► 6 and is indistinguishable from 5/2 in all simulated 
dimensions d > 6. But also the difference in the other ex- 
ponents between forest-fire and percolation values seems 
to vanish with increasing dimension. Of course a non in- 
teger upper critical dimensionality as d c = 11/2 would 
also be compatible with our results. 

C. The SOC forest— fire model on the Bethe lattice 

Often the critical exponents of a model in high dimen- 
sions are obtained not only by mean-field theory but also 
on the Bethe lattice. This might also be true for the SOC 
forest-fire model. The tree distribution on the Bethe lat- 
tice, however is not random and therefore different from 
mean-field theory, as the following consideration indi- 
cates: 
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A forest cluster of s trees on a Bethe lattice with con- 
figuration number z has 

s(z - 2) + 2 

empty neighbors, independently of its form. It is de- 
stroyed either when a tree grows at its edge (this hap- 
pens with probability p (s (z — 2) + 2) per time step), or 
when it is struck by lightning (with probability fs). A 
new cluster is generated each time a tree grows. In the 
steady state, the mean number of forest clusters is con- 
stant, i.e. [22] 

oo 

p (i - pi) = J2(f s + p ( 2 + s ( z ~ 2 ))M fi ) 

3 = 1 

OO 

= Pt(f + p(z-2)) + 2pJ2n(s). (4.5) 

s = l 

In the limit f/p — ► it follows from Eq. (4.5) 

oo 

J2n( S ) = (l- Pt (z-i))/2. (4.6) 

s = l 

If the stationary critical state were identical with mean- 
field theory of the forest-fire model, which is known to 
be identical with the mean-field theory of percolation 
[11,23], i.e. with percolation on a Bethe lattice, trees 
were randomly distributed with 

Pt = l/(z - 1) • 
With Eq. (4.6) this would lead to 

oo 
s = l 

i.e. the density of forest clusters would decrease to zero 
in the limit f/p — ► 0. Consequently, the exponent r were 
2 (see the comment after Eq. (3.8)). Both results are 
in contradiction to the well-known facts that a random 
tree distribution with density l/(z — 1) leads to a finite 
density of clusters and to r = 2.5 on the Bethe lattice. 
Therefore the SOC forest-fire model on the Bethe lattice 
has no random distribution of trees and is different from 
its mean-field-theory. But we cannot rule out the pos- 
sibility that the critical exponents nevertheless are the 
same as in mean-field theory. For the sandpile model 
on the Bethe lattice, it has been shown that the asymp- 
totic power laws are the same as in mean-field theory, al- 
though there exist short-range correlations between sites 
[24]. 

V. UNIVERSALITY OF THE CRITICAL 
EXPONENTS 

The critical behavior of a system usually depends only 
on properties as dimension and conservation laws, but 



not on microscopic details. We therefore expect that the 
critical exponents of the SOC forest-fire model are uni- 
versal under certain changes of the model rules. 

In order to check this assumption, we repeated the 
2D simulations for other lattice types. First, we chose a 
triangular lattice. This is equivalent to a square lattice 
with fire spreading to next-nearest-neighbors along one 
of the diagonals. Then, we investigated the model on a 
square lattice with next-nearest-neighbor interaction. 

The simulations of both variations of the model were 
done on a 4096 x 4096-lattice with f/p ranging from 
1/1000 to 1/8000. We compared the exponents r, p, p' , 
v, v' and a with the results given in Tab. I and found 
them to be exactly the same (see Figs. 7,8,9,10). 




1000 2000 4000 8000 

(f/p)" 1 

FIG. 7. The correlation length £ as function of (f/p) 1 . 
The slope yields the critical exponent v. (□ = square lattice, 
A = triangular lattice, * = next-nearest-neighbour interac- 
tion.) 




1000 2000 4000 8000 

(f/p)" 1 

FIG. 8. The lifetime of the largest fire T max as func- 
tion of (f/p) -1 - The slope yields the critical exponent 
v' . (□ = square lattice, A = triangular lattice, * = 
next-nearest-neighbour interaction.) 
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s 

FIG. 9. Fire lifetime T(s) for various lattice symmetries. 
The inverse slope yields the chemical dimension fi' . 




to 

FIG. 10. The temporal correlation function G(ui) for vari- 
ous lattice symmetries. The negative slope yields the critical 
exponent a. 

Looking at Eq. (3.1) one expects the forest density for 
the triangular lattice to drop below the square-lattice 
value of p c t 40.8%, and the forest density for next- 
nearest-neighbor interaction to drop even further, since 
the fire has now more possible paths to spread. The 
values obtained from our simulations are p'\ 34% for 
the triangular lattice and p c t 28% for next-nearest- 
neighbor interaction. 

Another modification of the model rules is obtained 
when trees are allowed to be immune against fire. A 
forest-fire model with immune trees, but without light- 
ning was investigated in [13]. Now we included the im- 
munity in the SOC forest-fire model by changing rule 2 
in the following way: 

• tree > burning tree with probability 1 — g n , if n 

nearest neighbors are burning. 

A detailed account on immunity in the SOC forest-fire 
model will be given elsewhere [25]. Here, we just state 
the main result: If one increases the immunity g from 
to some finite value below the critical value g c = 1/2, 



which is one minus the bond-percolation threshold, the 
critical exponents remain unchanged. In order to see the 
asymptotic critical behaviour in the simulation data how- 
ever, one has to go to very small values of f/p, when g is 
close to g c . For g = 1/2, the SOC fixed point disappears, 
and the size distribution of fires is identical to the size 
distribution of percolation clusters in uncorrelated bond 
percolation. 

VI. SUMMARY AND DISCUSSION 

In this paper, we have presented a scaling theory for 
the SOC forest-fire model, including also the dynamics 
of the system. The appropriate critical exponents were 
defined, and scaling relations between them were derived. 

The critical exponents in one dimension, which are 
known exactly [9], were rederived by simple arguments. 
They turned out to be just the classical ones proposed 
in [8]. In dimensions > 2, computer simulations then 
determined the values of the critical exponents and con- 
firmed the scaling relations. Our simulations have been 
performed closer to the critical point than any earlier 
simulation. The values of many exponents have been de- 
termined for the first time. Results which have already 
been obtained by other authors are mostly confirmed, 
some of them are corrected or improved. 

The simulations suggest that the critical exponents of 
the SOC forest-fire model in dimensions d > d c = 6 are 
given by its mean-field theory, which is identical with 
the mean-field theory of percolation. It still remains a 
challenge to give an analytic derivation of the upper crit- 
ical dimension. A short calculation showed that the SOC 
forest-fire model on the Bethe lattice is different from the 
mean-field theory of the model. 

Finally, we investigated the universality of the crit- 
ical properties by changing the lattice symmetry from 
"square" to "triangular" and by including next-nearest- 
neighbor interaction. Furthermore, the new parameter 
"immunity" was introduced into the SOC version of the 
forest-fire model. No violation of universality could be 
found. 

As already pointed out in earlier publications [13,14], 
there is a close relationship between the forest-fire model 
and excitable media which comprise phenomena so dif- 
ferent as spreading of deseases, oscillating chemical re- 
actions, propagation of electrical activity in neurons or 
heart muscles, and many more (For a review on excitable 
systems see e.g. [26,27]). These systems essentially have 
three states which are called quiescent (= tree), excited 
(= burning tree), and refractory (= empty site). Excita- 
tion spreads from one place to its neighbors if they are 
quiescent. After excitation, a refractory site needs some 
time to recover its quiescent state. In many of these sys- 
tems, spiral-waves have been observed. We expect that 
a SOC state can be found in some of these systems, if the 
appropriate range of parameter values is investigated, i.e. 
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if spontaneous excitation occurs rarely and if excitation 
spreads much faster than the system recovers from the 
refractory state. 
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TABLE I. Numerical results for the critical exponents in 1 to 8 dimensions (* = with logarithmic corrections, * = calculated 
from scaling relations). The exponents with index "perc" are those of percolation theory [19]. 



d 


1 


2 


3 


4 


5 


5 


7 


g 






1 6384 


448 


80 


32 


20 


1 2 


g 


T 


2 


2.14(3) 






2.45(3) 


2.50(3) 


2.50(3) 


2.50(3) 


"^perc 


2 


2 05 


2 18 


2 31 


2 41 


2 5 


2 5 


2 5 


A 


1 


1.15(3) 


1.30(6) 


1 ^ ( q ^ t 

1 .OO^O J 










1/6 


0* 


0.48(2) 


0.55(12) 












1/^perc 


1 


0.42 


0.56 


0.69 


0.85 


1 


1 


1 


r-2 

3 — T 





0.18 


0.30 


0.56 


0.82 


1 


1 


1 


Pi 


I 


0.4081(7) 


0.2190(6) 
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0.111(1) 


0.090(1) 


0.076(1) 


0.066(1) 


1 

2&-\ 


1 


0.333 
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1 96(1) 


2 51(3) 


3 


3 2(2) 








/iperc 


1 
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3.54 


4 


4 


4 




I* 


0.58 


0.52(3)^ 


0.53(3)^ 


0.57(7)t 








- 


1* 


0.43 


0.25 


< 0.02 t 


f 


f 


f 


f 


/ 

jl 


1 


1.89(3) 


2.04(10) 


2.02(10) 


1.98(10) 


1.94(10) 


1.92(10) 


1.93(10) 


jj,' 

Y 


1 


1.68 


1.89 






2 


2 


2 


y 


1* 


0.58 
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1.04(ll) t 


1.05(ll) t 
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~ I 

V 


1* 


0.44 


0.34 




o.i(i)t 


< 0.14 t 


< 0.08 1 " 


< 0.08 1 " 


z 


1 


1.04(2) f 


1.24(8) t 


1.49(10) t 


1.62(18) t 








b 


1 


1.27(7) f 


1.47(9) f 


1.73(10) t 


1.89(ll) t 


1.97(ll) t 


1.96(ll) t 


1.97(ll) t 


a 


2* 


1.72(5) 


2.15(5) 


2.00(5) 


2.01(5) 


1.95(10) 


1.96(10) 


1.96(10) 



11 



